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Abstract 

It is shown that the fluctuations of the jamming coverage upon Random Se- 
quential Adsorption ((Tgj), decay with the lattice size according to the power-law 
oq-j oc L -1 /^ 7 , with vj = 2D-d t > wnere D is the dimension of the substrate and 
(if is the fractal dimension of the set of sites belonging to the substrate where the 
RSA process actually takes place. This result is in excellent agreement with the 
figure recently reported by Vandewalle et al {Eur. Phys. J. B. 14, 407 (2000)), 
namely vj = 1.0 ± 0.1 for the RSA of needles with D = 2 and d{ = 2, that gives 
uj = 1. Furthermore, our prediction is in excellent agreement with different pre- 
vious numerical results. The derived relationships are also confirmed by means of 
extensive numerical simulations applied to the RSA of dimers on both stochastic 
and deterministic fractal substrates. 
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The irreversible deposition of particles on a surface involves two characteristic time 
scales: the time between depositions, and the diffusion time of the particles on the surface. 
For very strong interaction between particles and the substrate (chemical adsorption), dif- 
fusion becomes irrelevant and the Random Sequential Adsorption (RSA) model provides 
an excellent description of the underlying processes (for a review on RSA models see [1]). 
Under these conditions the system evolves rapidly toward far-from equilibrium conditions 
and the dynamics becomes essentially dominated by geometrical exclusion effects between 
particles. This kind of effects has been observed in numerous experiments [2]. 

The RSA of needles (or linear segments) on homogeneous, two-dimensional samples, 
has very recently attracted considerable interest [3, 4]. Particular attention has been 
drawn to the interplay between the jamming coverage and percolation [3, 4, 5]. The 
percolation problem has also attracted considerable attention in the field of statistical 
physics due to their relevance for the understanding of processes and phenomena in many 
other areas such as those occurring in disordered media, porous materials, systems of 
biological and ecological interest, etc. [6, 8, 9]. Therefore, a great progress in the field of 
the statistical physics of far-from equilibrium processes could be achieved by establishing 
links between RSA and percolation [3, 4, 5]. 

The percolation transition is related to the probability of occurrence of an infinite 
connectivity between randomly deposited objects, as a function of the fraction p of the 
substrate occupied by the objects. Close to the percolation threshold p c , the probability 
P to find a percolating cluster, on a finite sample of side L, is given by an error function 

[9] 

where a is the width of the transition region. It is well known that a vanishes in the 
thermodynamic limit according to [9] 

tJOcL - -, (2) 

where v is the exponent that governs the divergence of the correlation length as £ oc 
\p~Pc\~ u - 

Very recently it has been suggested that the jamming probability and the fluctuations 
of the jamming coverage may obey relationships similar to equations (1) and (2) [3], 
respectively. The aim of this note is to provide a qualitative derivation of equation (2) 
for the case of RSA on both homogeneous and deterministic fractal substrates. The 
predictions of the obtained equation will be compared with previously published data 
and further numerical tests will be performed. To accomplish these goals, the RSA of 
dimers on deterministic and stochastic fractals such as a Sierpinski Carpets (SC) [8] and 
the diffusion front [6, 8], has been studied. 

Let us first establish a link between the fluctuations of the number of deposited par- 
ticles (at the jamming state) on a subsystem of side L with those of a system of side 
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L, with Lq < L. Considering a fractal subsystem of side Lq, that has itself Qo minimal 
pattern blocks, and increasing the size of the subsystem n steps A times until reaching the 
size L, such as L(n) = \ n L , then Q will change to Q(n) = s n Q . Therefore, eliminating 
n, it follows that 



Q = Qox (^J (3) 

where df = log(s)/log(\) is the fractal dimension. 

Let Nq be the number of adsorbed particles in the starting subsystem of side Lq. For 
the system of side L the number of adsorbed particles N(L) is given by the sum 

N(L)=J2 N " W 
1=1 

where iVj are the number of adsorbed particles on each subsystem of side L that form 
the system of side L. Let (Jn(l ) be the fluctuations, in the starting L -subsystem, of N . 
If the correlation length associated to the random sequential adsorption (£r so ) is short 
compared with L (£r S(1 << L ), the random variables iVj will be statistically independent 
an so from Eq.(4) it follows 



s" 



'JV(L) 

i=l 



Furthermore, in the thermodynamic limit, the L -subsystems should have identical sta- 
tistical properties. Thus, their respective fluctuations will be the same. So, from the fact 
that s n = and Eqs.(3) and (5) one has 

= °k x L d \ (6) 
L o 

then the fluctuation of the density (9) in the system of size L can be obtained from Eq. (6) 
dividing by L 2D , so that 

i_ 

ere oc L -j, (7) 

where 



2D-df 

It should be stressed that Eqs. (7) and (8) are quite general relationships valid for 
substrate systems that are both homogeneous and deterministic fractal. Furthermore, the 
same relationships hold for the case of substrates globally-invariant under translations, 
such as random fractals, as it has been demonstrated elsewhere [11]. Also, the condition 
that the correlation length of the RSA process should be smaller than the system size 



3 



is usually valid for jammed states, where the correlation length is very short. It is also 
very interesting to notice that, using these relationships it may be possible to evaluate 
d { performing both RSA numerical simulations and actual experiments. Furthermore, 
existing numerical simulations performed in D = 2 dimensions with df — 2 are in ex- 
cellent agreement with equations (7) and (8) (notice that for these conditions it follows 
straightforwardly from equation (8) that vj — 1 exactly). In fact, for the jamming upon 
RSA of needles in two dimensions the value vj = 1.0 ±0.1 has been reported [3] and 
this figure is independent of the aspect ratio of the needles. Furthermore, early numerical 
results of Nakamura for the RSA of square blocks are also consistent with vj ~ 1 [10], 
while Kondrat et al. [4] have reported vj = 1.00 ± 0.05 for the RSA of segments on the 
square lattice. Since the obtained values for the exponent are independent (within error 
bars) of: i) the length of the segments (for all a = 1,2,.. ..,45) [4], ii) the aspect ratio of 
the needles [3] and iii) the size of the square blocks [10], it has been suggested that uj 
is a good candidate for an universal quantity of the jamming process [4]. Within this 
context, our finding shows that vj depends on the dimensionality of the substrate and 
the set where the RSA processes actually takes place. 

On homogeneous samples the jamming coverage (9) and its fluctuations (a$) can 
straightforwardly be obtained, since one has to deal with a single stochastic process. 
However, RSA on nonhomogeneous random substrates requires a careful treatment be- 
cause two correlated stochastic processes are now involved [11]. One can assume that the 
fluctuations due to the RSA process are given by an average over M independent samples: 



where a l e are the fluctuations measured using a single substrate sample but taken averages 
over independent RSA trials. It has been shown [11] that measuring a e with the aid of Eq. 
(9) one captures the physical behavior of the RSA process. In contrast, measuring the 
fluctuations of the average jamming coverage of different samples the physical behavior 
reflects the properties of the substrate [11]. 

In order to perform additional tests to the obtained analytical results, the RSA of 
dimers on both stochastic and deterministic fractals has been studied numerically. 

As example of an stochastic fractal, we have used a diffusion front. In order to generate 
the diffusion front, we considered the diffusion of particles at random, but with hard- 
core interactions, on a 2D square lattice of size L x L. There is a source of particles 
at the first row of the lattice y — 1,1 < x < L kept at concentration p(l,t) = 1. 
Also, at row y — L + 1,1 < x < L there is a well, p(L + l,t) = 0. So, there is a 
concentration gradient along the source-well direction, while along the perpendicular in- 
direction periodic boundary conditions are imposed. In the steady state the concentration 
gradient is constant, so one has 
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(9) 
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Vp{y)=L-\ (10) 

It is well known that the properties of the diffusion front [12, 13] are closely related to 
those of the incipient percolation cluster [6, 8, 9]. As the concentration p(y) of particles 
depends on the position, decreasing from the source to the well, one actually has a gradi- 
ent percolation system. The structure of the diffusion front is identical to the structure 
of the hull of the incipient percolation cluster [12]. Furthermore, the concentration of 
particles at the mean front position Uf is the same as the percolation threshold p c , so that 
P(Vf) = Pc [12]. The diffusion front is conveniently described by its average width o$ and 
the total number Nf of particles that constitute it. Using heuristic arguments it has been 
suggested that 

Nf 1 

-T-HVp(y/)r QJV where a N = — — (11) 

being v the critical exponent of the correlation length in the percolation problem [6, 8]; 
v = 4/3 in 2D, which gives a at = 3/7. So, from Eqs.(10) and (11) one has 

N f ~L d ?\ (12) 

with df F = 0:^ + 1 = 10/7^ 1.4286, and the diffusion front is a stochastic self-similar 
fractal [14]. 

RSA of dimers on diffusion fronts has been simulated using two rules: according to 
Rule I only adsorption events of dimers taking place on two nearest-neighbor (NN) sites, 
such us one of then belongs to the diffusion front and the remaining one is outside it, are 
considered. On the other hand, using Rule II one only allows the adsorption on NN sites 
of the diffusion front, disregarding adsorption trials on already occupied sites of the front 
and sites outside the fractal. 

RSA of dimers on deterministic fractals (Sierpinski Carpets [6, 8]) is also studied. The 
SC in D = 2 dimensions is generated by dividing a full square into \ D smaller squares 
of the same size. Out of these squares, k of them are chosen and removed. In the next 
iteration, the procedure is repeated by dividing each of the small squares left into A D 
smaller squares removing those k squares that are located at the same positions as in the 
first iteration. The resultant fractal dimensions are 

dt(s,\)=log(s)/log(\) (13) 

where s = \ D — k. In principle, this procedure has to be repeated again and again, 
however for the practical implementation in a computer only a finite number of iterations 
are actually performed [8, 7]. In a square lattice the smaller subdivision is actually a single 
site and the length is measured in site units. Furthermore there is a minimal pattern of A df 
sites. In the present work various generations of SC's of different size L, with periodical 
boundary conditions, have been employed. In all cases dimers are allow to adsorb only 
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Figure 1: Log-log plots of versus L for the case of RSA of dimers on the random fractal 
generated by diffusion fronts (DF). Results obtained using two different adsorption rules 
are shown. For details the adsorption rules see the text. 

on NN empty sites belonging to the fractal. For SC's with A = 3 and k = 1,2,3, 
as used in the simulations, the fractal dimensions are d^ = log(8)/log(3) ~ 1.8928, 
d fu = log(7)/log(3) « 1.7712 and rf flII = log(6)/log(3) « 1.6309, respectively. 

Figures 1 and 2 show log-log plots of a e versus L obtained upon RSA of dimers on 
diffusion fronts and Sierpinski Carpets, respectively. The obtained results, for these kind 
of fractals, are in excellent agreement with the prediction of Eq.(8) as follows from the 
comparison of evaluated and theoretical exponents listed in Table I. Further support to 
the theoretical prediction follows from additional results obtained using homogeneous 
samples, which are also listed in Table I. 

Summing up, it is shown that the exponent vj can be obtained as a function of the 
dimensionality D of the space and the fractal dimension df of the subset site where the 
RSA process actually takes place. Our main result vj = 2L f_ d , provides a solid ground 
to previous numerical data [3, 4, 10]. Furthermore, in this work, the validity of the 
proposed relationship is verified by means of extensive numerical simulations, using both 
homogeneous substrates as well as different fractals. 

Acknowledgments: This work was supported by CONICET, UNLP and ANPCyT 
(Argentina). 
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L 

Figure 2: Log-log plots of a e versus L for the case of RSA on Sierpinski Carpets obtained 
using different generating patterns as shown in the figure (black squares compose the 
fractal structure). 



Table 1: Examples of the application of Eq.(8) to different fractals as listed in the first 
column: SC = Sierpinski Carpet, DF = Diffusion front, HS2 Homogeneous Substrate in 
D = 2 dimensions. The 2nd column shows the exponents obtained fitting Eq.(8) to the 
simulation results while the 3rd one shows the estimations of dt obtained using — = 2D ~ d{ , 
The 4th column is a list of the exact values of df. Notice that for SC the labels a)-e) allows 
to identify the generating patterns, as shown in figure 2. 



Substrate 




dl 


df 


SC (a) 


1.051(4) 


1.898(8) 


ln(8)/ln(3) ~ 1.893 


SC (b) 


1.052(4) 


1.896(8) 


ln(8)/ln(3) ~ 1.893 


SC (c) 


1.115(2) 


1.770(4) 


ln(7)/ln(3) ~ 1.771 


SC (d) 


1.110(7) 


1.780(15) 


ln(7)/ln(3) ~ 1.771 


SC (e) 


1.16(2) 


1.68(4) 


ln(Q)/ln(3) ~ 1.631 


DF 


1.30(2) 


1.40(4) 


10/7 ~ 1.429 


HS2 (D = 2) 


1 




2 
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